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I. INTRODUCTION 



Particle creation in cosmological spacetimes was first discussed by Parker Q), Sexl and Urbantke (|], 
Zel'dovich and Starobinsky H in the late sixties. The basic mechanism can be understood as parametric 
amplification of vacuum fluctuations by an expanding universe . Particle creation in black hole spacetimes 
was first discovered by Hawking || (see also JjJ ) . A similar effect in a uniformly accelerated detector was 
discovered by Unruh g] and in a moving mirror by Davies and Fulling ^ . One special class of cosmological 
spacetimes which shows this characteristic thermal radiance is the de Sitter Universe, as shown by Gibbons 
and Hawking ]l0|] . One feature common to all these systems is that they all possess event horizons, and the 
conventional way to understand the thermal character of particle creation is by way of the periodicity on 
the propagator of quantum fields defined on the Euclidean section of the spacetime |ll],|l2| • 

One would not usually think of cosmological particle creation as thermal because in general such conditions 
(event horizon and periodicity) do not exist. However, several authors have shown that thermal radiance 
can arise from cosmological particle creation in spacetimes without an event horizon ||l3|-|l8|| . Each case has 
its particular reason for generating a thermal radiance, but there is not much discussion of the common 
ground for these cases of cosmological spacetimes. There also seems to be a gulf between our understanding 
of the mechanisms giving rise to thermal radiance in these two classes of spacetimes, i.e., spacetimes with 
and without event horizons^. 

In some of our earlier papers we have discussed thermal radiance in the class of spacetimes which possess 



event horizons (uniformly accelerated detectors 20 - 22 1 , moving mirrors and black holes 21 23j) using the 



viewpoint of exponential scaling |24|,|25| and the method of statistical field theory |26|,^7|. In a recent paper 
p3| we show how this method can be applied to spacetimes which possess event horizons only in some 
asymptotic limit, such as near-uniformly or finite-time accelerated detectors, and collapsing masses. In this 
paper, we study thermal particle creation from cosmological spacetimes with the aim of providing a common 
ground for cases where thermal radiance was reported before. Using the stochastic method, we show how to 
derive near-thermal radiance in spacetimes without an event horizon. 

The two primary examples we picked here for analyzing this issue are that of Parker and Berger [ p^pUl 
for an exponentially expanding universe, and that of slow-roll inflationary universes, in particular, the 
Brandenberger-Kahn model |18| . In Sec. 2 we examine several simple cosmological expansions which admit 
thermal particle creation and show their common ground, and their connection with thermal radiation in the 
(static coordinatized) de Sitter spacetime. We point out that all cases which report thermal radiation involve 
an exponential scale transformation |24|,^5|]. Thermal radiance observed in one vacuum can be understood 
as arising from the exponential scaling of vacuum fluctuations of the other vacuum pgfl . Sec. 3 applies this 
method to the Parker-Berger model to show how thermal radiance can be derived from exponentially-scaled 
vacuum fluctuations. Sec. 4 discusses particle creation from slow-roll inflationary universe and shows how 
near-thermal radiance can be derived in cases which depart from strictly exponential expansion. We end 
with a brief discussion in Sec. 5. 



II. THERMAL PARTICLE CREATION IN COSMOLOGICAL SPACETIMES: EXPONENTIAL 

SCALING 

Consider a spatially-flat (k=0) Robertson- Walker (KW) universe with metric 

ds 2 = dt 2 -a 2 Y,(dx 1 ) 2 , (2.1) 



1 Note that both of these effects are present in particle creation in a non-eternal black hole spacetime - over and 
above the thermal Hawking radiation for an eternal blackhole, there is also the contribution from backscattering of 
waves over a time-dependent classical effective potential, see, e.g. [^9) . 
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where t is cosmic time. A conformally-coupled massive (m) scalar field $ obeys the wave equation (e.g., 
©) 

[□ + m 2 + i?/6]$(i,x) =0, (2.2) 

where □ is the Laplace-Beltrami operator, and R — 6[a/a + (a/a) 2 ] is the curvature scalar. In a spatially- 
homogeneous space, the space and time parts of the wave function separate, with mode decomposition 
$>(t,x) — J^k 4>k{t)w k (x). For a spatially-flat RW universe Wk(x) — e lkx , and the conformally-related 
amplitude function Xk(v) = ufykif) of the k th mode obeys the wave equation in conformal time 77 = j dt/a: 

Xk (r))" + [k 2 + m 2 a 2 (r,)] Xk (v) = 0. (2.3) 

Call <f> l k n,ou (t, x) the modes with only positive frequency components at t- = —00 and t + = +00, respectively. 
They are related by the Bogolubov coefficients a k ,(3 k as follows: 

<P k n (t,x) = a k $l u \t,x) +/3 k ^*(t,x). (2.4) 

(For conformal fields it is convenient to use the conformally- related wave function X(rj,x) — a<f>(t,x). One 
can define the conformal vacua at rj± with ^ m > ou * i n terms of the positive frequency components.) The 
modulus of their ratio is useful for calculating the probability P n (k) of observing n particles in mode k at 
late times Q : 

P n {k) = \/3 k /a k \ 2n \a k \- 2 . (2.5) 

From this one can find the average number of particles (Nj:) created in mode k (in a comoving volume) at 
late times to be 

00 

n k = (N k ) = nP n {k) = \f3 k \ 2 . (2.6) 
n=0 



A. Bernard-Duncan Model 

In a model studied by Bernard and Duncan []l6| the scale factor a(rj) evolves like 

** Case 1 * * a 2 (rj) = A + Btanhprj (2.7) 

which tends to constant values a\ = A ± B at asymptotic times r\ — * ±00. Here p measures how fast the 
scale factor rises, and is the relevant parameter which enters in the temperature of thermal radiance. With 
this form for the scale function, a k ,(3 k have analytic forms in terms of products of gamma functions. One 
obtains 

\Pk/a k \ 2 — sinh 2 (7rw_/p)/smh 2 (7nj + / / 9) (2.8) 

where 

w± = i(w ou *±a; in ), < = ^ 2 + m 2 «l- (2-9) 

For cosmological models in which a(+oo) 3> a(— 00), the argument of sinh is very large (i.e. {tt/ p)u)± ^> 1). 
To a good approximation this has the form 

\l3 k /a k \ 2 = exp(-27ro; i „/p). (2.10) 

For high momentum modes, one can recognize the Planckian distribution with temperature given by 

k B T v = p/(2Tra+) (2.11) 

as detected by an observer (here in the conformal vacuum) at late times. 
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B. Parker-Berger Model 

This model is similar in spirit to the one proposed by Parker earlier E3], who considered a massless, 
minimally coupled scalar field in a Robertson- Walker universe with metric 

ds 2 = a 6 dr 2 - a 2 Y^{dx l ? (2.12) 

i 

where t is a time defined by dt = a 3 dr. The scale factor a is assumed to approach a constant at r — > 
±00, where in these asymptotic regions one can define a vacuum with respect to t time and construct the 
corresponding field theory. He considered the general class of functions for the scale factor (from Epstein 
and Eckart fl3lf ) 

**Case2** a 4 (r) = a\ + e CTr [(a 4 - a 4 )(e CTT + 1) + b] (e aT + 1)~ 2 (2.13) 

where a\,a2,b are adjustable parameters with a 2 > a\, and a is the rise function (similar to the p in the 
earlier case in terms of the conformal vacuum) . The modulus of the ratio of the Bogolubov coefficients is 
given by p) 



Oi- 



sin 2 ird + sinh 2 (irco- / a) 



(2.14) 



sin 2 ird + sinh 2 (7uj + /a) 

where d is a real number involving the constant b and 

u± = k(al±a 2 2 ). (2.15) 

For cosmological mod els in which a 2 3> ax, the argument of the sinh is very large, as in Case 1. Then, to a 
good approximation (2.14) has the form 

\Pk/a k \ 2 = exp(-4^fca 2 /a). (2.16) 

This form is independent of the adjustable parameters b, 0,2- From this one can show that the amount of 
particle creation is given by 

n k = [exp(fik) - I}' 1 (2.17) 

where [i = Ana\/a. Converting to the physical momentum at late times p = k/a 2 , one sees that this is a 
Planckian spectrum characteristic of thermal radiance with temperature 

k B T T = CT/(47ra?a 2 ). (2.18) 

How sensitively does the thermal character of particle creation depend on the scale factor? From physical 
considerations, the period when particle creation is significant is when the nonadiabaticity parameter f2 
satisfies J5) 

n = n'/n 2 >i. (2.19) 

Here f2 is the natural frequency (given by ha 2 in this case) and tt' — dVt/dr. Using this criterion, Parker 
argued that significant particle creation occurs during an early period when e T 1° (or in the first case, e v ^ p ) 
is small, whence a 4 has effectively the form 

* *Case 3 * * a 4 (r) = a\ + a^e T/ff (2.20) 

where ao is an adjustable parameter. This form of the scale factor was used by Berger |l3| ] for the calculation 
of particle creation in a Kasner universe, who also reported thermal radiation. Since particle creation vanishes 
at early and late times (as measured by the nonadiabaticity parameter), adiabatic vacua can be defined then 
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and one can construct WKB positive and negative frequency solutions for the calculation of the Bogolubov 
coefficients. Parker Q] showed explicitly th at th e modulus of their ratio has the same exponential form as 
that of the more complicated scale function (2.13), which, as we have seen, is what gives rise to the thermal 
character of the spectrum. Indeed he speculated that the exponential form in \@/a\ should hold for a general 
class of scale functions which possess the properties that: 1) they smoothly approach a constant at early 
time, 2) their values at late times are much larger than at initial times, and most importantly, that 3) they 
and their derivatives are continuous functions. The exponential factor contained in the scale functions in all 
three cases above at early times is thus responsible for the thermal property of particle creation, with the 
temperature proportional to the rise factor in the exponential function (a, or p in the first case). He also 
noted that this property is quite insensitive to the late time asymptotically static behavior of a(r) (could be 
the flattening behavior of a tanh function, or the rising behavior of an exp function). 



C. Common Features 



The three examples related above highlight an important common feature of thermal particle creation. 
That is, that a period of exponential expansion in the scale factor would give rise to thermal particle 
creation, i.e., an observer in the in vacuum before the expansion reports zero particles, while an observer in 
the out vacuum after the exponential expansion will report a thermal particle spectrum with temperature 
proportional to the rise factor in the exponential function in that particular time (e.g., conformal time rj in 
Bernard and Duncan's model, r in Parker's model, and cosmic time t in the de Sitter universe example below). 
The exponential scale factor and the relation between these two vacua are important to understanding particle 
creation in cosmological spacetimes on the same footing as that in the class of spacetimes with event horizons, 
including that of a uniform accelerated observer, a moving mirror, black holes and the de Sitter universe. 

To put the physics in a more general context, consider two vacua related by some transformation. Let 
us define the asymptotic in- vacuum as |0)t, the asymptotic out- vacuum as |0) otl t and the vacuum of an 
observer undergoing exponential expansion as |0) s . The in and out vacua are well defined because the 
scale factor approaches a constant at asymptotic past and future times, thus imparting the space with a 
Killing vector dt with respect to which one can define particle states in terms positive frequency modes. 
The s-vacuum is defined with respect to a different set of mode functions (like the Fulling- Rindler vacuum 
vis-a-vis the Minkowski vacuum for a uniformly- accelerated observer). The above examples calculate the 
particle creation between an in and out vacuum, but they also illustrate the important fact that the number 
of particles created is insensitive to the behavior of the scale function at late times - e.g., the result for the 
flattening tanh function which gives an asymptotically static universe is the same as that of the exponential 
function. Furthermore, the thermal nature of particle creation depends only on the initial stage of exponential 



expansion. In 1 28 2£j| these findings were used to connect the result of thermal radiance in these two classes of 
spacetimes. The assertion is that the more basic cause of thermality lies in the exponential scaling behavior 
rather than the existence of an event horizon [ p4|]25| ]. (The latter necessarily implies the former, but the 
converse is not always true) . The fine distinction between these two ways of understanding (kinematic versus 
geometric) thermal radiance will enable us to treat non-thermal cases in spacetimes which do not possess 
event horizons and to explore the stochastic nature of the Hawking- Unruh effect. 

At this point it is perhaps also useful for us to adopt this kinematic viewpoint to reexamine the cause of 
thermal radiance in a de Sitter universe. 



D. de Sitter Universe 

The de Sitter universe metric can be expressed in many different coordinates (see, e.g. f3(i|j ). There is the so- 
called closed (fc=l) Robertson- Walker (RW) coordinatization which covers the whole space, a(t) — cosh(Ht), 
the flat (fc=0) RW coordinatization which covers only half of de Sitter space with scale factor a(t) — e Ht , 
and the static coordinate shown below, to name just the commonly encountered ones. The vacuum states 
defined with respect to different coordinatization and normal mode decomposition have been studied by 
many authors p2| . In the static coordinate the metric is given by 
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d S 2 = [l-(gf) 2 ]d 2 - [i _ (Hg)2] . (2.21) 

Note that an event horizon exists at f = H' 1 for observers at f = 0, following the trajectory generated by 
the Killing vector d. This is similar in form to a Schwarzschild metric for a massive (M) object, restricted 
to 2 dimensions: 

, , 2AT , dr 2 , 

^ 2 = (!-— ) d ^- (Trwy- ( 2 - 22 ) 

for which Hawking first reported the famous black hole thermal radiance effect. 

Calculation of particle creation in the static de Sitter universe was carried out by Gibbons and Hawking 
JTcfl (GH) using the periodicity condition in the field propagator. Lapedes |33| gave derivations based on the 
use of Bogolubov transformations, and Brandenberger and Kahn treated the asymptotically de Sitter 
case. Let us analyze the relation between particle creation calculated in the Gibbons-Hawking vacuum 
(defined with respect to the 'static' de Sitter time ) |0) and that in the Robertson- Walker (k=0) vacuum 
(defined with respect to cosmic time t) \0)t- We will see that the |0) vacuum bears the same relation to |0)t 
vacuum as that between the exponential vacuum |0) s defined earlier and the asymptotic in- vacuum \0) t in 
the cosmological cases above. 



E. Exponential scaling: a kinematic viewpoint 

Starting from special relativity, assuming that two coordinate systems S = (t, r) and S = (, f ) (these 
are not the black hole coordinates) coincide at the origin, so the initial vacuum of S is the same as the 
RW vacuum, but the final vacuum in 5* is the GH vacuum. The two systems are related by the following 
conditions: 

(i) r = a(t)r (2.23) 
(it) a(t) = e m (2.24) 
(Hi) Hr = Har = dr = (3, (2.25) 

a(t) a/1-/? 2 

The meaning of these conditions is explained in p8| , which uses this example to illustrate the existence 
of exponential scaling in all cases which report thermal radiance. The two systems are related by a scale 
transformation (i), in this case, an exponential scaling (ii) such that an observer in S is seen as receding from 
S with a velocity of (3, with H the red-shift or Hubble parameter (in), and a Lorentz factor 7 (iv). Condition 
(iv) is called relativistic exponential transformation |2q|, w hich plays a central role in the understanding of 
the Hawking effect in terms of scaling concepts p^j25| , |34| |. With these correspondences, it is not difficult 
to see the analogy with the cosmological particle creation cases studied before. The initial (asymptotic in-) 
vacuum |0) f defined with respect to t time here (or tj,t time in the earlier examples) and the vacuum |0) 
defined in the exponentially receding system S bear the same relation. It is no surprise that the mod ulus o f 
the ratio between the Bogolubov coefficients have the same form characteristic of a thermal spectrum ( [2.10 ) 



but with p, a replaced by H . [One can find explicit calculation of thermal particle creation in Eq. (33) of 
p3| , using Bogolubov transformations replacing R by our H^ 1 .] The Hawking temperature for the de Sitter 
universe is given by 

k B T ds = (2.27) 

Once the relation between the de Sitter universe (which belongs to the class of spacetimes which show the 
Hawk ing-Unruh effec t), and that of some general cosmological spacetimes (with specific scale functions such 
as in (2.7, 2.13| , 2.20) is established, it is easy to generalize to the black hole and accelerated detector cases. 
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III. THERMAL RADIANCE IN THE PARKER-BERGER MODEL 



We now use the influence functional formalism (see e.g., p]f| ) to investigate particle creation in the Parker- 
Berger model [13 13. The line element is given by 



ds 2 = a 6 dT 2 -a 2 J2(dx 1 ) 2 



(3.1) 



where 



a 4 (r) = 1 + e 



(3.2) 



We consider the action of a massless, minimally coupled real scalar field 4>, which forms an environment 
acting upon a detector coupled to this field at some point in space. The field can be decomposed into 
a collection of oscillators of time-dependent frequency. Using the influence functional formalism, we can 
determine the effect of such an environment on the detector, which is also modelled by an oscillator. 
To do this we calculate the noise v and dissipation fi produced by the field. These are given by: 



( = v + ifi 



dk I(k,s,s')X(s)X*(s') 



(3.3) 



where / is the spectral density describing the system/environment interaction, and X is a sum of Bogolubov 
coefficients satisfying the classical equation of motion for the field oscillators. First we decompose the field 
into its modes; the Lagrangian density is 



C(x) 



-fJ 



In terms of normal modes the Lagrangian becomes 

1 



L(r)= £ ~ (O 2 - a^i&f 



k,a=± 



We see then that the bath can be described as a set of oscillators with mass and frequency 



m = 1 



4 1.2 



a 4 fc 



(3.4) 



(3.5) 



(3.6) 



Now as was already mentioned, X satisfies the classical equation of motion for an oscillator with the given 
parameters, and being a sum of Bogolubov coefficients, its initial values are predetermined. So we need to 
solve 



X"{T) + k 2 {l + e pT )X = 0, 

with initial conditions: 

X(r ) = l , X'( T0 ) = -ik. 
With a change of variables z = ln(2fc//?) + pr/2 we find the solutions in terms of the Bessel functions: 



X(t) = C1J2 



2k 



pr/2 



C2</_ 2>fc 



2k 



,PT/2 



(3.7) 
(3.8) 

(3.9) 
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To fix the constants ci,C2 consider that the initial time is To — * — oo; however the complex index Bessel 
functions oscillate infinitely often as their arguments approach zero, and so for now we leave To unspecified. 
In that case we can calculated c\ and C2; the final expression for X becomes, with 



2A- 
P 



f{r) = 



and Bessel indices labelled by v = 2ik/p: 



, „ irk , 2nk 1 

X(t) = — csch <^ ie pTo/2 

P P 



J v (J{r)) J-„(/(r)) 
J'M(ro)) J'-M(to)) 



J,A.f(r)) J-u(f(r)) 
Mf(r )) J-,(/(r )) 



(3.10) 



In the limiting case of t — * 00 (tq — > —00) we can use first order and asymptotic expressions for J and J' to 
write (which defines the phases a and 0) 



Mf(r)) 



I sinh 2tt k / p 
2-irk/p 



ia(k) 



-pr f) /2 / smh27rfc/p - m 

V 27Tk/ P 
2 ( . 7T 

— cos r(r) 

7TZ V w 2 4 



(3.11) 



In evaluating X(t)X*(t') we obtain various products of the Bessel functions with their derivatives (note: 
J* = J_ „); in particular we need 



/ 2ik\ (2ik 
[3 - a = arg T I H J - arg T I 

2ik 7r . 

= arg = — provided k f= 0. 

P 2 



(3.12) 



Also, when calculating the Bessel products, there arise sines and cosines with argument /(t) + f(r') = 
2k I p (^e pT / 2 + e pT / 2 ^j; when r — > 00 and we ultimately integrate over k, these terms won't contribute to the 
integral and so can be discarded. Changing to sum and difference variables defined by 



E = (r + r')/2 



A = t - t' 



we finally obtain 



dk I(k, t, t ) 



^fe^-e^) coth ^_ !Sln ^f 

£> V / O O \ 



cos 



2?rfc 
P 



2A- 



oP r/2 _ nr'/2 



We can now equate C with the standard form: 

/>oo 

C= / dkI eff (k,Y<) [C(/c,E)cosfcA-isinfcA] 
j 



(3.13) 



(3.14) 



(3.15) 



This is the form of ( which, with the function C replaced by coth would describe a thermal bath of static 
oscillators each in a coherent state. We will show that the unknown function C does indeed have the form 



2 To calculate these coefficients the wronskian of J2ik/ P an d J-nk/p is needed; note that there is a misprint in 
Gradshteyn and Ryzhik §8.474: the relevant quantity should be — — sini/7r. 
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of a coth, and can then deduce the temperature of the radiation seen by the detector. Here I e g(k, E) is the 
effective spectral density, also to be determined. We can always write C in this way since v is even in A while 
p is odd. By equating the real and imaginary parts of the two forms of £ and Fourier inverting, we obtain 



i r°° 

- dA cosfcA z/(E, A), 

K J-oo 

i r°° 

Ieff = / dA sinfcA /x(£,A). 

* J-oo 



(3.16) 



These expressions will be used throughout this paper to calculate C for the various cases of induced radiance 
that we consider. In order to use t hese w e need to calculate th e dis sipation and noise, p and v. 

We first evaluate p as given by (3.14); substituting it into (3.16) will then give us the effective spectral 
density I e ff{k, £). Define 

a = lLpa/2_ efU f/i\ = i e pS/2 smh P^ 

p V / p 4 

To proceed, we need to specify a form for the spectral density. This has been calculated in and in 3+1 
dimensions it is 



I(k,s,s') 



4^2' 



(3.18) 



where c is the coupling strength of the detector to the field. Then from ( 3.14 ) we have 

„2 



4tt 2 



- P 2/2 



dk k sin <rk 



An 2 



(3.19) 



where the last result follows from the discussion in |35[ |. Substituting this form for p into (3.16) gives the 
following result: 



c 2 k 

Ieff(k,X) = — e 



-3pE/2 



(3.20) 



Evaluating the noise kernel v is a more complicated affair. From ( 3.14 ) we write 



47T 2 



-pS/2 



dk k cos ok coth 



1-nk 



4tt 2 



d_ 

da 



(j 16 



1 P „_v.aAw' 



csch 



(3.21) 



where again the last integral has been calculated in |3J|]. Upon substituting this into (3.16) we obtain 

pS poo 



e pij r°° d 
C(fc,E) = — / dA cosfcA— Pill a) 
nk J_ 00 da 



167rfc 



dA cos fcA 



csc h 2 ^ - e" s csch 2 (e^' 2 sinh ^ 



(3.22) 



The first integral can be done by parts to get 
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dA cosfcA— P(l/a) = 
da 



dA d / / \ 

da — cosfcA— P(l/a) 
da da 



= -PV 

<y — oo 

47rfc 

= coth — 

P P 



dA d 
~a"d~A 

2nk 



dA 
da 



cos fcA 



(3.23) 



The second integral in (3.22) does not appear to be expressible in terms of known functions. Suppose we 
call it -B(fc, p, S), and write 



B(k,p,T,) = 2 dA cosfcA 
'o 



Then we have 



C(fc,/9,S) = e ps 



lcoth^ + /- e ^3(fc,p,£) 
p p Ibirk 



(3.24) 



(3.25) 



We need to examine the second term in the last brackets. The function B tends to zero for large k (by the 
Riemann-Lebesgue lemma), and attains a maximum at k = (since the cos term stops oscillating there). 
However, numerical work shows that this maximum value increases roughly with e pS which means that on 
first glance the second term in the brackets does not necessarily vanish at late times (S — > oo). So we need 
to examine the value of B at k = more closely, to see precisely how it changes with S. To this end we 



can consider -6(0, p, E) as a function of e pE and analyze its concavity; i.e. with x 



we need d 2 B/da 



Differentiating twice under the integral sign gives an integrand which is everywhere negative, and so we 
conclude that d 2 B/dx 2 < 0, which means that B as a function of x is everywhere concave down. But B 
increases with x, and thus Bjx = e~ pT, B — > as S — > oo. In that case the second term in the brackets gives 
no contribution in th e large time limit. 

Finally, from (3.15) we can write £ in a form which reveals the thermal nature of the detected radiation: 



C= / dfcI e# (fc,E) 
Jo 

The temperature of the radiation is then 



4e pE , 2nk , A . ' 
coth cos kA — i sin fcA 

P P 



(3.26) 



kuT 



4tt' 



(3.27) 



Inspection of (3.14) suggests that an alternative time variable can be chosen: 



The metric becomes 



P 



p 2 t 2 



(3.28) 



(3.29) 



with 



1 



2-I-2 



pH 



(3.30) 



Again following the previous formalism, we arrive at a description of the environment field in terms of 
oscillators, now with time dependent mass and frequency: 
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pt 

2 



4 1,2 



4a 4 fc 
P H 2 



(3.31) 



Solutions for X in this case are the same as before, and the calculations carry through in much the same 
way. With now S and A defined as mean and differences of t and t' we again arrive at thermal forms for 
the noise and dissipation: 



£ = e - P £/2 / dkl(k,s,s') 



cos fcA coth i sin fcA 

P 



(3.32) 



and the detected temperature is the same as ( |3 . 2 7| ) 



IV. INFLATIONARY UNIVERSE 



A. Eternal versus Slow-roll Inflation 



In this section we consider particle creation of a massless conformally coupled quantum scalar field at zero 
temperature in a spatially flat FRW universe undergoing a near-exponential (inflationary) expansion. The 
example of de Sitter space which corresponds to the exact exponential case has been treated in [^lj . Here we 
first solve for a general scale factor a(t) using a slightly different language from [^l). We then specialise to 
a spacetime (the Brandenberger-Kahn metric Jl8| ) which has initial de Sitter behavior but with scale factor 
tending toward a constant at late (cosmic) times. We can also define a parameter h which measures the 
departure from an exact exponential expansion. 

As before we first derive the noise and dissipation kernels by calculating X, the solution to the equation 
of motion of the field modes. The spatially- flat FRW metric is 



ds 2 = dt 2 -a 2 (t)J2(dx l Y 



(4.1) 



The Lagrangian density of the scalar field is 



C(x) 



(4.2) 



which leads to a Lagrangian in terms of the normal modes 



n 3 



k,a=± 



, it \ 2 CL „ _ f k a- \ , 

(€T + 1~&<& -[---) <T 

a \ or er 



(4.3) 



(We have added a surface term 1/a 3 d/dt (aa 2 q 2 ) to the Lagrangian. See |^T) for the rationale.) The classical 
equation of motion, and hence that of X, is 



X + 3-X 

a 



X = 



(4.4) 



with initial conditions 



X(U) = l , X'(t i ) = -ik-a'(t i ). 



(4.5) 



We find that 



X(t) = -e 
a 



(4.6) 
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where r\ is the u sua l conformal time. 

Now we use ( |3.3| ) to construct the influence kernel £, which contains the noise and dissipation kernels. 
From pl| the spectral density for the field is 



I(k,t,t') 



For the rest of this section we introduce for brevity: 



4tt 2 ' 



(4.7) 



x(E, A) = a(t)+a(t'), 
y(^,A) = V (t)-r 1 (t'). 



(4.8) 



Then we find from (3.3) 



( = v + ifi 



a(t)a(f) 4tt 2 J 



k e~ tky dk 



a(t)a(t') 4tt 2 



— P(l/y)+inS'(y) 



(4.9) 



As before we calculate the spectrum and temperature by Fourier transforming £ using ( 3.16 ). (Note that 
£, A are defined in terms of cosmic time t). We change the variable of integration from A to y, using 



dA 2a(t)a(t') 



(4.10) 



to write 



e 

e^k 
472 



S'(y)^dy 



(4.11) 



independently of the value of a(t). The temperature now follows from ( p.lq ): again we change to a y- 
integration by parts, remembering that A = y = 0. We obtain 



47T 3 



cos fcA dA d 
a(t)a(t') dy dy 



P(l/y) dy 



(4.12) 



which leads to 



C 



4 


/•OO 


d cos kA 


dA 


irk 


Jo 


dA x 


V 



(4.13) 



This equation is the central result of this section, in that it allows us to compute the spectrum corresponding 
to an arbitrary scale factor. For example, in the de Sitter case with a = e Ht we have 



= 2e H s cosh ■ 



HA 



2e~ H ^ HA 
y = H S ~2~ 



(4.14) 



which when substituted into (4.13) gives 



C = coth$. 

H 



(4.15) 



So for this case we can infer the temperature seen to be 



12 



(4.16) 



as was calculated with a slightly different approach in J2lj . 

As an aside, we note that from the above analysis for a general scale factor, the noise kernel is 



with dissipation 



f°° du 


d 


cos ku 


Jo y 


du 


ait) + a(t')_ 



e 2 5'(A) 

47T 



(4.17) 



(4.18) 



An often-used alternative to our principal part prescription is the introduction of a cutoff in the J °° expres- 
sions; unfortunately following this pro cedur e doesn't lead to tractable integrals even for the relatively simple 
de Sitter case. Note that in equation ( 4.13 ) for the temperature in the general case, we are essentially deal- 
ing w ith p roducts of a and r/, and it's therefore not surprising that for de Sitter expans ion, w here a oc 1/77, 
that ( 4.1S ) can be evaluated analytically. For other forms of a, even very simple ones, ( 1.15 ) becomes very 
complicated. 



B. Near-exponential expansion 

We now consider the case of a near de Sitter universe with a scale factor composed of the usual de Sitter 
one together with a factor that decays exponentially. We show that the spectrum seen is near-thermal 
tending toward thermal at late times. 

We start by considering the Hubble parameter to have a constant value (characterising de Sitter space) 
plus an exponentially decaying term: 

Hit) = Ho (1 + ae-P Hot ) , (4.19) 
and from this our aim is to calculate C using (|4.13| ). The scale factor results from integrating H, and is 



a(t) = exp ( Hot 



-(3H t 



We dehnc the parameter h which measures the departure from exact exponential expansion to be 

-afie"^ 



hit) = H{t) 



Hit) 2 

as fit — > 00, and as we might expect it is exponentially decaying at late times. 

To proceed, we indicate the de Sitter quantities by a subscript zero as well as writing 



(4.20) 



(4.21) 



E = HoT, 



A = HoA. 



Then immediately we have from (4.14) 



xo = 2e cosh ■ 



yo 



2e- s . , A 

sinn — . 

H 2 



(4.22) 



(4.23) 



We wish to perturb these by using the new scale factor. Suppose we define perturbations fx, / 2 by writing 

x = xJl + f 1 {%'k)) , y = y (l + f 2 (%A)). (4.24) 
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We first have 



x = a(t) + a(f ) = e Hot -f e ^ Hot + e ^'~f^" ot ' . 



which in the late time limit can be approximated by 

2a 



which yields f± 



Next we write 



. _ -/JE C ° Sh 2^ 

J 1 ~ 6 



cosh 



A 



v = v(t) - v(t') = 



dt 
, ajt) 

and by making the same late time approximation as for x we get 

ri 



y 



-Hot 



1 



a 



-PHat 



This leads to 



2ae -(i+/3)s _ (i + / 3)A 

= Vn H ; ; Sinn . 

y ° (3(1 + 0)H O 2 



ae -/3S sinh (lM)A 

J2 — 



+ sinh A 



Note that at late times fa, fa tend to zero. In that case to calculate C we write ( 4.13 ) in the form 



C 



4 f°° d 
irk J dA 



tkA 



Xq 



(1 - h) 



1-/2 

yo 



dA, 



and so write, to first order in fa, fa: 



C 



irk 



dA d cos fcA 
yo dA x 



dA 
itk J y 



fa 



d cos fcA d fa cos fcA 



dA xq 



dA xq 



=cothTrk/H =AC, the perturbation 

Evaluating AC is lengthy but straightforward so we merely write the answer in integral form: 



AC = 



27rfc/? 



dA 



-2k sinh 



(1+/3)A 



sin fcA 



1 sinh 



(1+/3)A 



;fcA 



Ho(l+P) sinh 2 # C osh f l + P cosh 2 A sinh f 



2fc cosh (1 ~ fflA sin fcA (1-/3) sinh (1 ~ fflA cosfcA 2 cosh (1 ~ fflA cos fcA 



H ° cosh 2 f sinh f 



cosh 2 — sinh A 



cosh 3 f 



The important point is that the factor e " ensures that this perturbation to the thermal spectrum 
exponentially at late times. 
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C. Brandenberger-Kahn model 



We are now in a position to derive the function C(k, E) for the Brandenberger-Kahn model. In this case, 

o(t)=e ^(l-e-/>) (4 34) 

with Ho, a constants. As t tends toward zero and infinity, a(t) tends toward e Ht and e 2H ^ a respectively. 
The Hubble expansion function is 



H(t) = - = H a e- at ' 2 (4.35) 
a 

and the parameter which measures the departure from exact exponential expansion h(t) is 



Hit) 



at/2 



2H 



2H Q 



+ 0(a 2 t 2 ) 



(4.36) 



We assume that \at\ <C 1. Eqn ( 4.1 3| ) is much too difficult to evaluate analytically here, but we can get some 
insight by calculating it as a first order correction in h to the de Sitter case. 

At this point, we also mention an alternative perturbation of de Sitter space, given by the scale factor 



fjA P H(t)dt 

a(t) = e-io y ' 



(4.37) 



which describes a solution of the vacuum Einstein equations with a time-dependent cosmological constant 
A(t) = 3H 2 (t). One may expand H(t) in a power series about t = 0. Defining h as in (4.36), this form of 
perturbation turns out to be identical to the Brandenberger-Kahn model to first order in h. We have, to 
first order, 



H(t) = H Q + H 2 ht, 



(4.38) 



the correspondence between h and a being given by ( 4.36 ). We will therefore calculate the detector response 
for the Brandenberger-Kahn model only, keeping in mind its correspondence with the model mentioned 
above. 

Again define f\ , f% (now with h included) such that 



x = x + hf 1 (A) , y = y + hf 2 (A). 
The corrections are then written as 



(4.39) 



A(A) 



/ 2 (A) = - — 

-no 



AM A A 

— cosh h EA sinh — 

4/2 2 



~„ A 2 ~ \ A 
E H — + 2E + 2 sinh — 



1 A cosh — 
2 



After some computation we obtain the spectrum to be 



C(fc,E) = (l + ZiTi) coth 



a form which shows its approximately thermal nature, with 



nk 



Ti = - 



E^ 
sinh # 

no 



cosh 



r + (E + l) 



1 — tanh 



irk\ 2 



2ku 
Hi 



Hq J ft Jo sinh 



^ du 



(4.40) 



(4.41) 



(4.42) 
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As a function of k/Ho the integral looks much like tan 1 , tending to tt/2 as k/Ho — > oo. 
In the low frequency limit the departure from a thermal spectrum is, to 0(k 2 ): 



hTi ~ ft 



I - ,v + 2/8) (g 



(4.43) 



Note that we stipulated that \hT\ ~ |ai| <C 1, so that /iTi remains small as time passes. In the high frequency 
limit the departure is given by 



hT, -» -2ftS e-'^a ^ _ ^ 
which again remains small, and is especially close to zero for high frequencies. 



(4.44) 



V. DISCUSSION 



This paper continues the theme of our previous papers on the stochastic approach to particle creation 
pl| ^3],^9) with focus on two main points: 1) A unified approach to treat thermal particle creation from 
both spacetimes with and without event horizons based on the interpretation proposed by one of us [ p4)|25| ] 
that the thermal radiance can be viewed as quantum noises of the field amplified by an exponential scale 
transformation in these systems (in specific vacuum states) fl28|| . In contradistinction to viewing these 
as global, geometric effects, this viewpoint emphasizes the kinematic effect of scaling on the vacuum. 2)An 
approximation scheme to show that near-thermal radiation is emitted from systems in spacetimes undergoing 
near-exponential expansion. We wish to demonstrate the relative ease in constructing perturbation theory 
using the stochastic approach. 

The emphasis of the statistical field theory is on how quantum and thermal fluctuations of the matter fields 
are affected by different kinematic or dynamic conditions. For particle creation in spacetimes with event 
horizons, such as for an accelerated observer and black holes, this method derives the Hawking and Unruh 
effect p(j| , f2lf from the viewpoint of amplification of quantum noise and exploits the fluctuation- dissipation 
relation which measures the balance between fluctuations in the detector and dissipation in the field [ p2[ . 
For spacetimes without event horizons, such as that in near-uniformly accelerated detectors or collapsing 
masses and wide classes of cosmological models, some studied here, one can describe them with a single 
parameter measuring the deviation from uniformity (acceleration) or stationarity (expansion) which enters 
in the near-thermal spectrum of particle creation in all these systems. The fact that we can understand all 
thermal radiation generating processes in these two apparently distinct classes of (cosmological and black 
hole) spacetimes ^§^| and be able to calculate near-thermal radiance in this and earlier papers testifies to 
the conceptual unity and methodological capability of this approach. 
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